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Abstract 



We present a new class of matrix models which are manifestly symmetric under the T- 
duality transformation of the target space. The models may serve as a nonperturbative 
regularization for the T-duality symmetry in continuum string theory. In particular, it now 
becomes possible to extract winding modes explicitly in terms of extended matrix variables. 
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I. INTRODUCTION 



Matrix models can be regarded as a nonperturbative regularization of string theories. 
Recently, it became clear that various duality symmetries play vital roles in understanding 
the nonperturbative properties of string theories. Unfortunately, however, it does not seem 
that the standard matrix models are particularly suited for studying the duality symme- 
tries. Even the T-duality which is valid to each order of genus expansion in continuum 
string theories is difficult to implement manifestly in matrix models, at least in their present 
interpretation. For example, as discussed in our previous paper [|lj , the Ising model described 
by a standard two-matrix model does not preserve T-duality symmetry, whose validity is 
naively expected from the Kramers- Wannier duality on a square lattice, once we take into 
account the higher genus effects. The reason for this failure is that there is no symmetry 
between global winding and momentum modes in such models. Although it is not difficult 
to define partition functions as sums over random surfaces preserving T-duality [], it is quite 
nontrivial to write down corresponding matrix models preserving T-duality symmetry man- 
ifestly. Since at present matrix-model approaches seem to be the only known tractable way 
towards dynamical formulation of nonperturbative string theories, it is worthwhile to develop 
methods of treating duality symmetries exactly using matrix models from various possible 
standpoints. This is particularly so in view of a recent trend concerning new possibilities of 
using matrix models in seeking for fundamental theories of strings including D-branes. 

The purpose of the present paper is to present a new class of extended matrix models 
which have exact and manifest T-duality symmetry under the usual random-surface inter- 
pretation of matrix models. Although the models in general are not exactly soluble using 
presently available techniques, we believe that the existence of such models is interesting 
by itself and will be useful for future investigations. To maintain the T-duality symmetry 
exactly, we need two key ingredients that have to be combined in a suitable way. One is that 
the global, as well as local, degrees of freedom corresponding to the winding and momentum 
modes must appear symmetrically for arbitrary topology of surfaces. The other is that the 
random triangulation of surfaces must be invariant under the duality transformation. For 
the first property, we adopt a variant of Z K spin systems (K = an integer) as the target 
space, generalizing the discussion given in the Appendix A of our previous paper jL]. For 
the second, we can use the result of reference || which is, to the best of our knowledge, the 
only known matrix model with exact self-duality. In this reference, the dual transformation 
for the random triangulation (namely, the measure of pure 2D gravity) was discussed, but 
the first problem was not addressed. In fact, a naive generalization of the method used for 
the one-matrix model does not work for systems coupled with matter. We will clarify how 
to combine self-dual matter systems with the method of 13]. 

The plan of the paper is as follows. In the next section, we briefly discuss a class of 
Zk spin systems on a fixed surface, which we call Zp t Q models (K = PQ) and have man- 
ifest symmetry between "momentum" and "winding" modes analogously with the toroidal 
compactification of continuum string theory. In section III, we present our extended matrix 
models by coupling the Zpq models to 2D gravity with the self-dual measure, and exhibit 



See, e.g., ^ in which the case of c = 1 has been discussed. 
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their exact T-duality symmetry. In section IV, some simple cases will be discussed for the 
purpose of concrete illustration. We will then discuss some relevant issues including the 
c = 1 limit. 



II. Z PtQ MODEL ON A FIXED SURFACE 

As is well known, Zk spin systems in two dimensions 0] are self-dual under the Kramers- 
Wannier dual transformation. There is, however, a subtlety when we do not neglect global 
degrees of freedom. To make the system completely self-dual, we have to impose appropriate 
constraints in order to suppress the local vortex excitations which tend to violate the duality 
symmetry. In the case of the lattice c = 1 matter system, this was discussed in and the 
Appendix A of [Jj . Here we present its generalization to Zk systems. For clarity of notation, 
we assume the surface to be a fixed square lattice and denote the spin sites by x and links 
by x, [L where \i is the direction of the link from a site x. The difference operation between 
the nearest neighbor sites x and x + fi where fi is a unit lattice vector along the fi direction 
is denoted by A M . 

For a given K, we introduce two integers P, Q satisfying K = PQ. Spin variables 
n x — 0, 1, . . . , P — 1 (mod P) are then assumed to live in the space Zp. In addition to the 
spin variables, we introduce link variables m 1](1 = 0, 1, . . . , Q — 1 (mod Q) living in the space 
Zq. Then the partition function of our Zpg model is defined as 

Z = (R E )IF (Q) ( A /^- A u m x ^l[B(A^n x -Pm x ^. (2.1) 

x,[i n x ,m x ^ x x,)i 

The Boltzmann factor B(x) is assumed to be periodic under the translation x — > x + K, 
and the 5-function constraint is understood modulo Q. Note that after the summation over 
the link variables, the Boltzmann factor can be regarded as a function on Zp with respect 
to the spin variables. Obviously, if the topology of the fixed surface is sphere, the constraint 
leads to 

m x ^ = A^m x , m x G Z Q . (2.2) 

Therefore by redefining the spin variable by n x — > n' x = n x — Pm x G Zpq = Zk, the system 
is reduced to an ordinary Zk spin system. If the fixed surface is not the sphere, we have 
instead of flO) 



m Xifl = A^m, + m x ^, (m x , m x ^) e Z Q , (2.3) 

where m Xjfl is a global vector field which cannot be reduced to the difference of the site 
variables and is associated with a nontrivial homology cycle of the surface. The degree of 
freedom represented by m x ^ is the analogue of the winding modes in the toroidal compact- 
ification of continuum string theory. 

Now let us perform the dual transformation on the partition function ( |2.1| ). We introduce 
two auxiliary fields n y G Zq and ij) y il G Zk on dual sites y and dual links, respectively, and 
rewrite (|2.1|), apart from a numerical proportional factor, as 



(n e n e )ne i[ ^^ A ^ + ^ e '^ (A, ' n ~ pm '' )i n^(^,M)- (2.4) 

x,[i n x ,m XilJ , y,^ h y ,ip ytll x y,[i 
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Here we suppressed the subscripts for the sites on the exponential to avoid unnecessary 
notational complexity. The Boltzmann factor B is the Fourier transform of the original 
Boltzmann factor: 

B{a)= ' Ee'fr*B(6). (2.5) 

2^b=0 B \ b ) 6=0 



Reality of the Boltzmann factor leads B(b) = B(K — b). We assume the normalization 
condition B(0) = B(0) = 1. Then solving the constraint coming from the summation over 
m M , the general solution for ip^ is 

^Pfi, = A M n - Qrh^, rh^ e Z P . (2.6) 

After substituting this result and taking the summation over the original spin variables n, 
we have the constraint, 

e^A>, = (modP). (2.7) 

Thus the partition function after dual transformation takes, apart from a numerical overall 
normalization factor, the same form as the original one (|2.1| ) with the interchange P <-> 
Q, B B 

Z = ([[ E )U 5(P) ( A ^- A^)n^( A A-K^ ( 2 -8) 

y,n h y ,m y ^ y y,n 

The global modes represented by fh y ^ are interpreted as the analogue of the momentum 
mode of toroidally compactified strings, and hence the dual transformation interchanges the 
winding and momentum modes, precisely as required for T-duality symmetry. In particular, 
the systems with P = Q are self-dual on surfaces of arbitrary genus. 

Note that, although the models with a given K = PQ are equivalent to the usual Zk 
spin models on the sphere, they are in general different on higher-genus surfaces, because 
of the different appearance of the global momentum and winding modes, depending on the 
choice of P and Q. In a sense, the Zp t Q model amounts to compactifying the Zk target 
space by a subgroup Zp of the isometry group Zk- For example, the standard Ising model 
corresponds to P = 2, Q = 1, Z 2i ± model. Hence its dual is Z 12 . Therefore the Ising model 
is not exactly self-dual on higher genus, as discussed in |]]. In this class of models, the 
simplest exactly self-dual model is Z 2t 2, which is identical with the (Z 4 ) Ashkin- Teller model 
on the sphere. Finally, we note that taking the limit P = Q — > oo appropriately gives a 
lattice version of toroidally compactified strings discussed in the Appendix A of A brief 
discussion on this limit will be given in section IV. 



III. SELF-DUAL MATRIX MODELS: GENERAL THEORY 

We now show how to construct the matrix model corresponding to the Zp^Q model. In 
the standard method of coupling the Z^ spin system to 2D gravity, we introduce K different 
Hermitian N x N matrices M a , a = 1, 2, ... ,K and assume the action 
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S = NTy 



1 K K 

^ o,6=l o=l 



(3.1) 



In the Feynman graph expansion, the potential V(M a ) represents a spin site which is on the 
center of a discretized surface element. The propagator C~^, on the other hand, corresponds 
to the Boltzmann factor assigned to a link connecting the nearest neighbor spin sites with 
spin variables a and b. The Zk symmetry requires that the kinetic operator (and the 
propagator) satisfies translation invariance modulo K with respect to the indices a, b, 

C ab = C(a -b) =C(a-b± K), C^ 1 = D(a -b) = D(a-b± K). 

In our terminology, this construction only represents the Zk,i model. To extend the 
construction to general Zpq matrix models, we first introduce a set of P different Hermitian 
QN x QN matrices, denoted as M a , a = 1,2, ... ,P. The index a ranging from 1 to P 
corresponds to the Zp spin variables as above, while the additional Q x Q matrix elements 
of each M a are supposed to be associated with the Zq link variables. To prevent unnecessary 
complication in notations, we will always suppress the original U(N) indices in the following 
and only indicate the additional Q x Q matrix indices by e {1, 2, ... , Q}. The choice 

of an appropriate propagator will then enable us to construct the Zpq matrix model. The 
correct choice is, suppressing the U(N) indices, as follows: 

D Hm = TfJl B ( a - b - Pm){L m UL^)kv (3-2) 

iV m=0 

where B is nothing but the Boltzmann factor used in the previous section, and L is an 
arbitrary Q x Q unitary matrix satisfying the conditions 

L Q = J, L- 1 = L\ (3.3) 



tr Q Z/ = if % ^ (modulo Q). (3.4) 

Here the trace notation tiq means taking the trace only with respect to the indices 
The notation Tr will be used for the total trace operation including both the U(N) and the 
U(Q) indices. On the other hand, the trace operation with respect only to the ordinary 
U(N) indices will be denoted by tr^y. A formula for the inverse of the propagator Q3.2|) will 
be given later. 

The potential function corresponding to the spin sites are assumed to be the following 
special form using the result of [|J 

p 

NTiJ2H^-gM a ). (3.5) 

a=l 

This is necessary to ensure that the measure for the random triangulation is invariant under 
the T-duality transformation even before taking the (double) scaling limit. 

The model is invariant under the group U(N) x U(Q),M a -> UM a U~ l , U G U{N) x 
U(Q), provided the matrix L is transformed by the U(Q) part of U. It is easy to check that 
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the conditions ( |3.3|) , (|3.4j ) precisely impose the constraint corresponding to the 8^> function 
in ( |2.1| ) for general triangulation of arbitrary surfaces, at each elementary closed circle of 
links (i.e., plaquette). Hence, the matrix L is arbitrary under the U(N) x U(Q) invariant 
conditions ( |3.3| ), (|3.4|). Also it is not difficult to prove explicitly that this model restricted to 
the sphere approximation is equivalent to the standard Zk matrix model as briefly discussed 
in section IV. 

Now that we have given the definitions of the matrix models, let us next proceed to 
rewriting of the models in a manifestly T-duality symmetric form. We will extend the 
method given in || for the case of pure gravity to our Z P q model. For this purpose, it is 
convenient to go to a particular representation for the matrix L using the U(Q) symmetry, 
namely the diagonalized representation given by 



Lij = Sije'v^. (3.6) 



Thus the propagator is now 



D a i'u = ^ £ B(a-b- Pm)5 u 5 kj e^- j)m . (3.7) 



-i 



m=0 



Here, by putting additional Kronecker S's, 5 ac , 6bd, we duplicated the indices a and b in order 
to make the appearance of the i, k (g Zq) and a, b (g Zp) indices symmetric. We note that 
in this representation the propagator has manifest Zq periodicity with respect to the indices 
while the periodicity with respect to the indices a,b,... is not manifest. This is 
expected from the situation of the fixed-surface model of the previous section, since the Zp 
periodicity appears only after summing over the link variables. In the matrix model, the 
latter operation appears only for Feynman amplitudes after taking the trace operation. 

We first introduce PN auxiliary vector fields ipf which transform as a complex vector 
under the group U (N) x U (Q). Here as before we have suppressed the U(N) vector ( "color" ) 
indices and also the "flavor" U(N) vector indices, while both of the matter indices a(G U(P) 
vector) and z(g U(Q) vector) are explicitly indicated. Thus the auxiliary fields i[)f can 
actually be treated as NQ x NP matrix fields which transform as 

< -> (U c ^V f )t, 

where U c G U(N) x U(Q) and Vj G U(N) x U(P). Using these auxiliary fields, rewrite the 
potential term as 

e -mr N J2 P a=lH i- gMa ) = J dW t e -ELxE2^^ tB a-^^. (3.8) 

Here we used the notation ip* for the complex conjugate of the ^'s treating them as N x iV 
(color x flavor) matrices. We can then perform the integration over the matrix M a . Using 
the propagator ( p.7|) , the partition function takes the form 

Q P 1 Q P 



Z = N p ,q{B) J dipdtp^ exp 



££tr^r + ^ 2 E E tr^lX-M 



2 

i=l a=l i,j,k,l=l a,b,c,d=l 



(3.9) 
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The normalization constant Np,q(B) originates from the Gaussian integration, 



Mp,q{B) = det'/'D, 



where D is the matrix of the propagator 

The T-dual transformation amounts to rewriting the partition function in terms of the 
dual Boltzmann factor B(a), defined by ( f2~5[) . By directly substituting (|2.5| ) to the propa- 
gator (|3T7j), we find 

D$£t = AT ^K Q if, (h , E B(i - J - Qmje'K-X-^W^^^^, (3.10) 

A little examination of this expression shows that the partition function in terms of the B 
takes the same general form as the original one apart from the normalization factor, 



Z = M P>Q (B) J d4>d$ exp 



EX>^M* + ^ 2 E E teNM'BS&tiW 

a=l i=l a,b,c,d=l i,j,k,l=l 

(3.11) 



after interchanging the Zp and Zq indices by making the following redefinitions of the 
auxiliary vectors, propagator and the coupling constant, respectively, as 



iA\ (3-12) 



= ^E B(l ~ J ~ Qrh)e^ d -^5 ab 5 cd 5 H 5^ (3.13) 



m=0 



9 2 Q = g 2 Y^ B{a) 

z^m p 



9 2 ^~9 2 = " ir% jr. = ^= " ■ (3.14) 



Here we have interchanged the color and flavor indices by performing a cyclic permutation for 
the auxiliary fields ?/>'s in the quartic term of the action 0, corresponding to the replacement, 

Di%1(B)^e^-^Wt%(B), 

of which the phase factor is absorbed by ( |3.12| ). We note that in this dual form, the Zp 
periodicity instead of the Zq periodicity is now manifest. This comes about due to the 
redefinition (|3.12|) of the auxiliary fields: The field ipf is supposed to be manifestly Zq 



i — i 

periodic with respect to the index i, while it acquires a phase e Q under the translation 
a — > a + P. This property is interchanged after the above redefinition from ip to <fi. It should 
also be remarked that the auxiliary fields ip an d 4> carry both the matter (Zp) and dual- 
matter (Zq) indices in addition to the color and flavor U(N) indices. This is an interesting 



2 Note that we have used the same trace notation tr^r for U(N) flavor indices. 
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aspect of the present formalism, which is something not appeared in the familiar formulation 
of T-duality in string theory. 

The above result is precisely the dual transformed structure of the Zp^q models, apart 
from the redefinition of the coupling constant which is connected to the overall normalization 
factor suppressed in the previous discussion for the fixed surface. It is evident that by 
reversing the route from the M a (a = 1, 2, . . . , P) representation to the ipf representation, 
we can construct the dual transformed matrix model, now with Q Hermitian NP x NP 
matrices Mj (i = 1,2, ... ,Q) from the 0„ representation. Thus we have established that 
our matrix models have desired T-duality properties. Formally, the T-duality symmetry 
corresponds to the following identity for the partition function, 



l P , Q (B,g) _ Z Q , P (B, : 



Mp,q{B) MqAB) 



(3.15) 



Note that the normalization constants Mp t Q and jVq,p whose ratio is easily calculable using 
( p.lOj ) are independent of the coupling constant g. As a consequence, the normalization 
constants can be neglected in the scaling limit. 

In order to write down the actions of the matrix models explicitly in terms of the Boltz- 
mann factor B, we need a general formula for the inverse of the propagator ( |3.2|) . This is 
easily obtained by applying Zp Fourier transformation appropriately. The result is 



Cab / j-) — l\ab 

ij,kl — \ U )ij,kl 



1 P- 1 1 



NSaS^e-^-^ £ t e <g(«-ft)™ (3 . 16) 
P m=0 B(m;i - k) 



where 



p-i Q-i 

■ 2n 



B{m; i - k) = Y, E B ( a ~ P/) e -^ (a -" )(l - A;) e-^ ma . (3.17) 

a=0 (=0 

Corresponding to the Zq periodicity of the propagator, this form of the kinetic term is Zq 
periodic with respect to the indices while under the Z P translation it is periodic 

only up to a phase e~ l ~$^~ k \ Of course, by inverting this formula, we could have started 
from the general form of the kinetic term having the same periodicity and expressed the 
Boltzmann factor in terms of the kinetic term. 

Finally, we briefly touch upon the question of observables in our models. Since we have 
to preserve the U(N) x U(Q) symmetry, the set of the most general invariants consists of 
the traces of arbitrary polynomials consisting of the matrices M a and £ where £ is the 
U(N) x U(Q) matrix acting as the identity in U(N) and as L in U(Q): 

A ab ... (n, /,...; m, p, . . .) = ^Tr(M a ™£ m M^/7 • • •). 

If the total power of the matrix £ is q (mod Q), A a j,... (n,l, . . . ;m,p, . . .) represents a loop 
state with "winding number q" . In the standard models, it has been impossible to explicitly 
extract the winding modes in terms of the matrix variables. Our construction thus indicates 
how to remedy this deficiency of the usual matrix models by extending the U(N) symmetry 
to U(N) x U{Q). 



S 



IV. EXAMPLES AND DISCUSSIONS 



The purpose of this section is to first present a few simple special cases of our extended 
models in more concrete forms as an illustration of the general theory, and then further 
discuss some of the important aspects of the models, in particular, the reduction of the 
degrees of freedom and the c = 1 limit. 



A. Examples 

We give concrete formulas for three cases, Z 12 , Z 1>3 and Z 2)2 - 

(i) Duality between Z\^ 2 and Z 2) \ matrix models 

This corresponds to the case treated in |5],[l| with cubic potential. The Z\$ model has 
only link variables. The dual version, on the other hand, is the Z 2 \ model which is nothing 
but the ordinary two-matrix model with spin variables, but without any link variables. Let 
us confirm this by an explicit computation. The propagator of the Z\ 2 matrix model is 
written as 

{M^Mm) = -(5 u 6 kj + e~ 2(i L a L kj ), (4.1) 



where (5 is an inverse temperature of the Z 12 spin system. L has the properties in Eqs. (ft.3|) 
and ( |3.4| ) with Q = 2. Correspondingly, the action is 

S M = N - ^ 4 ^ Tr(M 2 - e~ 2/3 £M£M) + iVTr ln(l - gM). (4.2) 

Introducing the auxiliary variables ipi and performing the M-integral, we have the action 

^ = N^N^i - N^g 2 tv N J (£ + e " 2/3 ( E 4U^f \ ■ ( 4 -3) 

i=l [ i=l i,j=l ) 

For the purpose of direct transition to the dual-matrix representation in the simple example 
treated here, it is convenient to introduce two auxiliary matrices U and V, which are Her- 
mitian and have U(N) flavor indices, instead of faithfully following the procedure adopted 
in the general theory. Then the action can be rewritten as 

S m = Nti N 4^ + nUi n {U 2 + e^V 2 ) - Ng%T N (Y,^\ + E ^^4). (4.4) 

i i,j 

Integrating ip and changing the basis as 

M 1 = U + V, M 2 = U - V, (4.5) 

we get the dual action 

1 1 



S M = N— _-t rjv 



.V/ 2 • .V/.7 - 2e - ; .\7,.\'/ 2 +Nj2^NH 1 ~9 M a), (4.6) 



a=l 
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where we rescaled M a as M a — > y/l + e _2/3 M a corresponding to the redefinition ( |3.14j ). The 
dual temperature f3 is defined by the Fourier transformation (|2.5|) as 

e~ 2 ^ = tanh (3, 
which is of course the famous Kramers- Wannier relation, 
(ii) Z 13 matrix model 

We repeat a similar calculation for the matrix model. The propagator and the action 
have the following forms: 



(MijMki) = —[SuS kj + e- 3 ^L u (L 2 ) kj + e-^(L 2 ) a L kj ], 



(4.7) 



l + 2e"^2 



M — 



1 + ea 



-CMC 2 M 



+ iVTrln(l -gM). 



(4.8) 



Through the same steps as before, we obtain the action of the variables ^'s. Then, 
introducing a Hermitian matrix U and a complex matrix X, we rewrite as 

3 i 

Suxt = iV^tr^M + N-tr N (U 2 + 2ef /3 XtX) 



-NgJ2tr N (AU4)-Ng £ tr^X^jL^ + ^}(L%^X). 



i=l 



After the ^-integral and the replacement 

M x = U + X + X\ M 2 = U + uj 2 X + ujX\ M 3 = U + ujX + u 2 X\ 
with ijj = e l ~ , we obtain the action of the Z 3 1 matrix model 



S M = N- 



cothp 1 



-tr 



N 



l + 2e"^2 

3 

+Nj2^ N \n(l-gM a ), 



Ym 2 — 



{M X M 2 + M 2 M 3 + M 3 Mi 



0=1 



(4.9) 



(4.10) 



(4.11) 



where the rescaling M a — > y 1 + 2e 2^M a was done, and the dual temperature is given by 



e 2 ' 



C2 



^ + 2 



C2 



(iii) Z 2 $ matrix model 

This is the simplest example of self-dual models. If we adopt the standard Boltzmann 
factor for the model, 



^n x - 2m X]tl ) = e^ > 



/3[cos^-(A t ,n x ~2m x ,^-l] 
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having a self-dual structure with 



e^ = coth^, (4.12) 



the propagator and the corresponding action are given, respectively, as 

((MO^MOh) = ((M 2 ) ij (M 2 ) kl ) = ^nSkj+e-^LaLkj), 

((M 1 ) ii (M 2 ) w ) = le-PfaSu + L u L kj ), (4.13) 

e 2/3 i ( 2 ^| 

S m = N - . ^ 2 -^ Et( M -) 2 + e- 2f3 CM a CM a ] - 2e-?M 1 M 2 - 2e" /3 £M 1 £M 2 
4smh (3 2 [ a=1 J 

2 

+AT^Trln[l - gM a \. (4.14) 

a=l 

Following the steps of the general theory, we can easily confirm that the model is self-dual 
with respect to ( |4.12j ) and g = g 1+ ^ ' ■ 



B. Reduction of degrees of freedom 

From the discussion of section II, it is clear that the Zpq matrix models must be equiva- 
lent to the standard Z& matrix models when restricted to the sphere approximation. Appar- 
ently, however, the Zpq matrix models have more degrees of freedom than the models 
with K = PQ. If counted as Af x JV Hermitian matrices, we have PQ 2 matrices instead of 
PQ. Let us therefore examine how the reduction of degrees of freedom occurs. It will turn 
out that there is a sort of local gauge symmetry in Feynman diagrams which accounts for 
the superfluous degrees of freedom. 

For simplicity, we here consider a concrete example of the Zp 2 (Q = 2) models. Once 
we understand this special case, extension to the general case will be straightforward. For 
this purpose, it is more advantageous to use the real representation of the matrix L given 
by 

Lij = S%h, (4-15) 

than the diagonal representation used in the general theory in section III. The nonzero 
components of the propagator then satisfy the following symmetry properties with respect 
to the Z 2 indices (i, j type) 

-Dll.ll = D 22 ,22 = -Dl2,21 = -021,12, 

-Dll,22 — -022,11 — D 12 i 2 = -021,21) (4.16) 

where we have suppressed the Zp indices (a, b type) since the identities are valid for arbitrary 
a, b. All other components of the propagator vanish. In particular, there is no coupling 
between the diagonal and off-diagonal matrix elements with respect to the Z 2 indices. This 
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shows that the links assigned to the off-diagonal matrix elements must always form Z 2 
closed loops. Furthermore, an amplitude with the off-diagonal elements is identical with 
a corresponding diagram in which all the off-diagonal matrix elements are replaced by the 
corresponding diagonal elements as indicated in (|4.16|) . Note here that the vertices behave 
as identities with respect the Z 2 indices and hence do not have any contribution violating 
the symmetry between the diagrams with diagonal and off-diagonal elements. In the case 
of sphere topology, the correspondence between the diagrams is established by making the 
interchange 1 <-> 2 globally for all Z 2 indices inside the domain enclosed by the closed curve 
of off-diagonal links. This replaces the off-diagonal links along the closed curve to links with 
diagonal matrix elements, without changing the amplitude. We can start this procedure 
from the smallest domains and go to larger ones successively to eliminate all off-diagonal 
links in this way. Thus the Feynman diagram contributions have a kind of gauge degeneracy 
2 L which is determined by the number L of independent Z 2 closed loops of off-diagonal links 
on the surface. Using the fact that L coincides with the number of dual sites, it is easy to 
see that the free energy becomes essentially identical with the model with only the diagonal 
matrix elements (M a )a, (i = 1, 2) under the rescaling of N and the coupling constant as 
N — > 2N,g n — > 2 _ 2 _ g n where g n is the coupling constant for the term ^TrM™ in the 
potential. Thus the matrix spin degree of freedom is reduced to 2P which is the same as 
for the standard Z 2 p. Obviously, the present argument cannot be extended to higher genus, 
since the above correspondence between diagonal and off-diagonal amplitudes cannot be 
established when the closed loop of off-diagonal links wraps around a non-trivial homology 
cycle. 

Extension to the general case is straightforward. The links with off-diagonal elements 
form Zq closed loops and the propagator satisfies the identities D^^i = D i+m j +n ^ k+n i +m . 
Repeating the above arguments, we arrive at the rescaling N — > QN, g n — > Q~*~ g n after 
the system is reduced to matrices with diagonal (with respect to the Zq indices) elements. 
The same result can also be obtained by a simple counting of the Feynman diagrams by 
focusing on their dependence with respect to N, Q and the coupling constant. 



C. c = 1 limit 



We next turn to the problem of the c = 1 limit of our models which corresponds to 
taking the limit P = Q — > 00. This problem is important if one regards the present matrix 
models as a nonperturbative regularization of the T-duality symmetry of string theory. 

First of all, take the Zq^q model assuming the most general form of the Boltzmann factor 
H as 



B(a) 



exp 




2-kol5 



cos ■ 



Q 2 



1 



where -KVs are constants. In the limit Q — > 00, we obtain the standard Gaussian Boltzmann 



factor as follows: when Kg is large with 
obtain 



j3s being fixed, we can expand the cosine to 
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R 2 2vr 
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a - 27rQn M ) 



(4.17) 



In the last expression, we left one of /Vs with 5 = 0(1) nonzero and redefined y~ = /3,5<5 2 to 
make B(a) fit to the ordinary Gaussian Boltzmann factor of the c = 1 model. 
For taking the c = 1 limit on a fixed lattice, we make the following replacement 



2nn 271^ 



2tt 



dX. 



Then the partition function of the Zq,q model becomes, apart from a numerical proportional 
factor, 



r2n 

Z = J dx S^A^exp 



N u ez 
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47T 



(A M X - 2^ 



(4.18) 



where we rewrite m M + Qn^ by N^. Here and below, we suppress the subscripts for sites 
and correspondingly the product symbols for sites for notational simplicity. In the dual 
representation, corresponding to n £ Zq and ip £ Z^ in section II, it is natural to take the 
limit as, 

2vrn ~ 2tt _ r** ~ 2ir _ /■«> 

Similarly to the section II, the summation over leads to ^ = A^X — 2-nm^, fh^ £ Z, 
and the summation over imposes the constraint on m M . Thus the partition function after 
the dual transformation takes the same form as (|4.18 ), 
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Z (x. \ dX X ^A M m„exp 



rhu,eZ 



^(A U X - 27cm u ) 2 



(4.19) 



with the correspondence R «-> 1/R. As has been discussed in |2|] and the Appendix A of 
U it is also straightforward to show self-duality by a direct dual transformation ( [4.18j ) <-> 
( |1.1S ) in the continuous target space. 

Let us next consider the Q — > 00 limit of the matrix model corresponding to the ^q,q 
model. Since the model is manifestly self-dual for an arbitrary integer Q, all we need to 
check is whether the dual transformation laws of physical parameters are well-defined in 
the limit. There are two such parameters, namely, the compactification radius R and the 
coupling constant g. First, as for R, it follows from ( f4.18| ) and fl4.19p that R is transformed 
as R <-> 1/R, since the transformation law of the Boltzmann factor is of course identical 
with the case of the fixed lattice. The transformation of the coupling constant is given by 
( p,14| ) which in the present limit leads to 



~2 9 

9 =R- 



(4.20) 



Thus both are well-defined in the c = 1 limit. 
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In particular, nonsingular transformation property of the coupling constant would allow 
a well-defined double scaling limit in the sense of world sheet. From the relation ( 4.20Q , we 
expect that the critical point function of R would behave as 



g c = f{R)R l '\ 

where f(R) is a function invariant under the dual transformation R «-> 1/R. Unfortunately, 
however, self-duality alone is not powerful enough for obtaining exact results in the scaling 
limit. For example, it is possible that the function f(R) actually vanishes in the c = 1 limit 
due to the degeneracy, as suggested from the counting argument of the previous subsection. 
In that case, the interplay between the c = 1 limit and the scaling limit is rather subtle, 
since taking the formal zero-coupling limit g — > directly in the action is not meaningful and 
we have to examine the physical correlators as in the usual situation in general lattice field 
theories. In any case, however, our discussion is sufficient to show that the present models 
can serve as a nonperturbative regularization of the T-duality symmetry in the continuum 
critical string theories. In particular, it is guaranteed that the properly defined continuum 
partition function is duality symmetric, as has been explicitly computed Q in the standard 
c = 1 model of finite radius. 

D. Conclusion 

We have established a new class of extended matrix models which have manifest sym- 
metry under the T-duality transformation. We can think of several possible extensions and 
applications of our work. Among others, we would like to mention the following. 

1. Construction of manifestly T-duality symmetric string field theories: Using the macro- 
scopic loop variables including the L matrix, we can try to derive stochastic Hamilto- 
nians following our previous works |J[T|]. 

2. Extension to open strings: In particular, it is interesting to formulate some toy models 
for Dirichlet branes using the present method. It might be useful for studying the 
nonperturbative dynamics of the D-branes. 

3. Inclusion of fermions: In particular, extension to supersymmetric models. 

4. Application to critical strings by extending the target space to higher dimensional 
spaces: For this, however, we have to deepen our understanding of the matrix-model 
formulation of critical strings. 

In all of these possibilities, it would be more or less crucial to develop some powerful methods 
for treating the scaling limit in order to perform useful nonperturbative analysis on the basis 
of our models. We hope to return to some of the above issues in future works. 
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